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ABSTRACT 

It has been suggested that the peculiar properties of the luminous Type Ic supernova SN 1998bw and its low- 
energy gamma-ray burst GRB 980425 may be understood if they originated in a standard gamma-ray burst explo- 
sion viewed far from the axis of the relativistic jet. In this scenario, strong radio emission is predicted from the 
jet on a timescale 1 to 10 years after the explosion as it decelerates and spreads into our line of sight. To test this 
hypothesis we have carried out late-time radio observations of SN 1998bw at t = 5.6 years, yielding upper limits 
which are consistent with the continued fading of the supernova. We find these limits to be consistent with an 
off-axis jet only if the progenitor mass loss rate is M < 4 x 10~ 7 M yr" 1 (for a wind velocity v w = 1000 km s" 1 ) 
or the fraction of the shock energy in magnetic fields is £b < 10~ 3 . These values are low relative to those inferred 
for cosmological GRBs. We combine the SN 1998bw measurements with existing observations for a sample of 
15 local Type Ibc supernovae to estimate that at most 6% produce collimated, relativistic outflows. 

Subject headings: supernovae: individual (SN 1998bw) - gamma rays: bursts 



1. INTRODUCTION 

GRB 980425 was the first gamma-ray burst (GRB) to be 
identified with a supernova (SN), SN 1998bw (Galama et al. 
1998; Pian et al. 2000). However, this association was not 
universally accepted since the spatial and temporal coincidence 
of these events could only be approximated to ±8' and ±2 
d, respectively. Recently, the detection of spectroscopic fea- 
tures in the light curve of GRB 030329, similar to those seen 
in SN 1998bw (Stanek et al. 2003; Hjorth et al. 2003), has 
strengthened the hypothesis of the SN 1998bw/GRB 980425 
association. These results have given a powerful impetus to 
the "collapsar" model (MacFadyen & Woosley 1999; Zhang, 
Woosley & MacFadyen 2003) in which a Wolf-Rayet progeni- 
tor undergoes core collapse, producing a rapidly rotating black 
hole surrounded by an accretion disk which injects energy into 
the system and thus acts as a "central engine". The energy ex- 
tracted from this system supports a quasi-spherical Type Ibc SN 
explosion and drives collimated jets through the stellar rotation 
axis which produce the prompt gamma-ray and afterglow emis- 
sion (see review by Zhang & Meszaros 2003). 

Despite this progress, some unresolved issues remain. Fore- 
most among these is understanding the connection between lo- 
cal SNe Ibc and cosmological GRBs. Estimates of the fraction 
of Type Ibc SNe that produce a GRB range from 10~ 5 (Perna, 
Sari, & Frail 2003), to 0.5% (Frail et al. 2001), and even ap- 
proaching 100% (Lamb et al. 2003). This uncertainty in the 
relative rates of Ibc SNe and GRBs is due to different assump- 
tions about the geometry and energetics of GRBs. 

As the only local SN observed in association with a GRB, 
SN 1998bw is the key to our understanding. At a distance 
of 38 Mpc, this Type Ic SN was unusually luminous at opti- 
cal and radio wavelengths (Kulkarni et al. 1998; Galama et 
al. 1998). The broad absorption lines seen in early spectra 
of SN 1998bw implied expansion velocities of > 30,000kms _1 



and (isotropic) kinetic energies of ~ 3 x 10 52 erg (Iwamoto et 
al. 1998; Woosley, Eastman & Schmidt 1999; Tan, Matzner & 
McKee 2001). Likewise, the bright, early-peaked radio emis- 
sion implied a significant amount of energy (~ 10 50 erg) cou- 
pled to mildly (T ss 2) relativistic ejecta (Kulkarni et al. 1998). 
Detailed modeling by Li & Chevalier (1999) confirmed this re- 
sult and also showed the need for a second energy injection 
episode, indicating the presence of a central engine, similar to 
the model inferred for GRBs. 

However, in contrast to the large energies inferred from the 
optical and radio emission from SN 1998bw, GRB 980425 was 
a sub-energetic gamma-ray burst. The prompt emission had an 
(isotropic) energy release of only 8 x 10 47 erg (Pian et al. 2000) 
- 4 to 6 orders of magnitude below typical GRBs (Frail et al. 
2001). 

One popular explanation posits that SN 1998bw/GRB980425 
was a typical GRB viewed away from the jet axis (Nakamura 
1998; Cen 1998; Eichler & Levinson 1999; Woosley, East- 
man & Schmidt 1999; Salmonson 2001; Yamazaki, Yonetoku 
& Nakamura 2003). This hypothesis can been described by 
two different scenarios based on the ratio of the off-axis an- 
gle, # oa , to the opening angle of the GRB jet, 6j, In Case 1, 
# oa ~ 39 j so the jet emission is detected at early time. Here, 
the 7-rays originate from the edge of the relativistic jet, caus- 
ing the inferred (isotropic) energy to be suppressed (Granot et 
al. 2002). In Case 2, 6 0a is large (>> 10 Oj) so the jet emission 
is only detectable at late time (t <~ 1 — 10 years) when the jet 
has reached spherical symmetry. In this case, the prompt 7-ray 
emission may be due to Compton scattering of photons into our 
line-of-sight (Waxman 2004a). 

Granot et al. (2002) have investigated Case 1, finding the op- 
tical and 7-ray emission to be consistent with an off-axis angle 
of 9 oa « 4° (for T w 100). This scenario, however, may have 
difficulty explaining the X-ray and radio evolution. Waxman 
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(2004a) has recently investigated Case 2, predicting a late-time 
rise in the observed luminosity which is most easily detectable 
at radio frequencies (Paczynski 2001; Totani 2003; Granot & 
Loeb 2003). 

The possibility of Case 2 has given rise to unification mod- 
els which imply some fraction of local Ibc SNe can also be 
described as GRBs viewed off-axis (Waxman 2004a; Lamb et 
al. 2003). While the off-axis jet model provides a convenient 
framework in which to unite the GRB and SN phenomena, con- 
firmation of the model requires observational evidence for a 
GRB jet within a local SN. In this paper we carry out late-time 
radio observations of SN 1998bw in an effort to detect the puta- 
tive off-axis jet from GRB 980425. We combine this measure- 
ment with existing observations for a sample of nearby Type 
Ibc SNe to place constraints on the parameters of the off-axis 
jet and the fraction of local Ibc SNe which could be associ- 
ated with relativistic jets similar to those seen in cosmological 
GRBs. 

2. OBSERVATIONS 

We observed SN 1998bw with the Australia Telescope Com- 
pact Array (ATCA) 4 beginning at 2003 December 4 21:00 UT 
and ending on 2003 December 5 8:30 UT. Dual frequency ob- 
servations were made at 1384 and 2368 MHz, using a band- 
width of 128 MHz and two orthogonal linear polarizations for 
each frequency band. The final synthesized beams were 8 x 14" 
and 5x8", at 1384 and 2368 MHz, respectively. 

No emission is detected toward SN 1998bw at either fre- 
quency but there is a weak source (209 yuJy) detected at 2368 
MHz 1.7" west and 3.1" north of the SN. This source is too 
bright to be due to thermal emission from the complex of HII re- 
gions in this area. It is possible that this is the radio counterpart 
of the ultraluminous X-ray transient (S lb) recently identified by 
Kouveliotou et al. (2004) but it would be 60-times more lumi- 
nous than the only other known detection (Kaaret et al. 2003). 
The peak flux density at the position of the SN was 62 ± 40 
and 149 ± 50 /zJy for 1384 and 2368 MHz, respectively. These 
new upper limits, taken 2049.19 d after GRB 980425 exploded 
(t ~5.6 yrs), lie just below the extrapolation of a power-law de- 
cay from 50 to 790 days (Figure 1). 

3. SN 1998BW: CONSTRAINING THE OFF-AXIS JET MODEL 

Waxman (2004a) has made predictions for the evolution of 
late-time radio emission from a typical GRB jet viewed far 
from the jet axis. As the off-axis jet sweeps up material and 
decelerates, it eventually undergoes a dynamical transition to 
sub-relativistic expansion (Frail, Waxman & Kulkarni 2000) 
and spreads sideways into our line-of-sight. This occurs on a 
timescale f oa « E$\/(Ms/vw.3) years. Here £51 is the energy of 
the jet (normalized to 10 51 erg), v w is the velocity of the stellar 
wind, and Ms /v Wj 3 is the mass loss rate of the progenitor star 
normalized as (M/10 -5 M yr _1 )(v w /10 3 km/s) and equivalent 
to A* = 1 in the notation commonly used for GRBs (Chevalier, 
Li & Fransson 2003). On this same timescale the radio emis- 
sion from the SN, which is produced by a spherical shock in- 
teracting with the circumstellar medium, will fade significantly 
with a power-law decay index of a m 1 .3 - 1 .6 (Li & Chevalier 
1999). 

Particle acceleration within the jet produces bright syn- 



chrotron emission which causes a rapid rise of the radio lu- 
minosity as the jet comes into our line-of-sight. Making the 
typical assumption that the energy is partitioned between the 
electrons and the magnetic field, and that these fractions (e e and 
eg) are constant throughout the evolution of the jet, the spectral 
luminosity is given by Waxman (2004a) as: 



U, = 



1.7 x 10 30 (3e c )(3e B ) 3/4 x 
(M. 5 /v W 3p 4 / v \-i ~ 
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On timescales t < f oa the spectral luminosity of the jet rises 
steeply as f +3 9 (Granot et al. 2002) and during t > f oa it de- 
cays as f~ 1,5 , assuming an electron spectral index p = 2. As 
indicated by Eqn 1, an off-axis GRB jet is most easily detected 
at radio frequencies. 

Late-time radio observations of SN 1998bw can therefore be 
used to impose constraints on the parameter space for an as- 
sociated off-axis GRB jet. In modeling the evolution of the 
putative jet emission, we consider four model parameters, £51, 
M-5/v w ,3, e e and e B , for which we adopt values comparable to 
those inferred for for cosmological GRBs. The kinetic energy 
of the jet and the energy fraction in electrons are fairly well 
constrained by broadband observations of GRB afterglows, so 
we assume typical values for these parameters of £51 = 1 and 
e e = 0.1 (Frail et al. 2001; Panaitescu & Kumar 2001; Yost 
et al. 2003). The energy fraction in magnetic fields is signif- 
icantly less constrained, however, showing a typical range of 
values between e# = 0.002-0.25 (Yost et al. 2003). This range 
imposes a significant uncertainty in the luminosity of an off- 
axis jet, through Eqn 1 . 

The model is also strongly dependent on the mass loss rate of 
the progenitor - generally believed to be a massive Wolf-Rayet 
(WR) star. Mass loss rates inferred from broadband model- 
ing of GRBs are typically M_ 5 /v Wj 3 « 0.1- 1 (Chevalier, Li & 
Fransson 2003; Panaitescu & Kumar 2002; Yost et al. 2003), 
yet local WR stars have observed mass loss rates significantly 
larger, on the order of Ms w 0.6-9.5 (Cappa, Goss & van 
der Hucht 2003). The uncertainty associated with the progen- 
itor mass loss rate maps to a large range in the predicted peak 
time and luminosity of an off-axis jet. Thus the uncertainty in 
Ms I v w ,3 and e# produces a degeneracy in the model which we 
are able to constrain using late-time radio observations of local 
Ibc SNe. 

3.1. The e B - Ms / vy 3 Degeneracy 

Figure 1 shows the predicted radio luminosity evolution of an 
off-axis GRB jet based on Eqn 1 and using values of Ms / vy 3 = 
0.01,0.1, 1 and 10. This plot demonstrates how the uncertainty 
in €b produces a family of lightcurves for each Ms /vV3 value. 

The full set of constraints is derived through investigation 
of the two-dimensional parameter space of e B -Ms/v w ^. Fig- 
ure 2 shows how the predicted spectral luminosity of the pu- 
tative jet at time, t , maps to a single contour in this parameter 
space. By comparing the observed luminosity of SN 1998bw 
at each epoch with the jet luminosity predicted for that time, 
we exclude the region of parameter space rightward of the cor- 
responding contour since this region produces a jet which is 
brighter than the observed SN 1998bw emission at that epoch. 
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The region of reasonable parameter space for the putative jet 
is bracketed by M_5 / v w ,3 = 0.04- 6 as derived from modeling of 
SN 1998bw radio light curves (Li & Chevalier 1999; Waxman 
2004b) 5 , and by e B = 0.002-0.25 as inferred for most cosmo- 
logical GRBs. The summed constraints (from t = 1 1 to 2049 
days) rule out the majority of this region. 

4. LOCAL TYPE IBC SUPERNOVAE: FURTHER CONSTRAINTS 

The search for an off-axis jet can similarly be carried out to- 
ward other local (di < 100 Mpc) Type Ibc SNe for which there 
are late-time (t ~ 1-20 year) radio observations. Eight type Ibc 
SNe were taken from Berger et al. (2003a): seven upper lim- 
its (SN2001B, SN2001ci, SN2001ef, SN2001ej, SN2001is, 
SN2002J, SN2002bm), and a detection (SN 2002ap) at t w 
0.5 years. We have supplemented this sample with Very 
Large Array 6 archival observations of SN 1983N, SN 1984L, 
SN 1985F, SN 1990B, SN 19941, and SN 1997X taken at t w 
16.7, 2.9, 7.1, 0.8, 8.0 and 2.4 years, respectively. In addi- 
tion, we include the recent detection of SN2001em at t k, 2.4 
years (Stockdale et al. 2004) - the only SN within this sample 
for which there are no early time radio observations. The SNe 
data are plotted in Figure 1. With the exception of SN 2001 em, 
all of the late-time observations are significantly fainter than 
the M-s/v Wi 3 = 1 off-axis jet prediction and nine SNe constrain 
M- 5 /v w , 3 <0.1. 

Using the same method applied in §3, we derive constraints 
for these 15 SNe. The resulting contours are plotted with 
SN 1998bw in Figure 3. For SNe with later observations, the 
constraints on e B values improve, as demonstrated by the con- 
tours for SN 19941 and SN 1983N. The more robust constraints 
are provided by the faintest luminosity limits given by the near- 
est type Ic supernovae, SN 19941, SN2002ap, SN 1983N and 
SN 1985F, all at <; 8 Mpc. Additional constraints are derived 
from the SNe for which we detect late-time emission, SN 19941 
and SN 2002ap 7 , since they show evidence that the radio lumi- 
nosity is decaying with respect to earlier lightcurve measure- 
ments (see Berger, Kulkarni & Chevalier 2002) and are thus 
inconsistent with rising off-axis jet emission. 

The region of reasonable parameter space for an off-axis 
jet is roughly bracketed by M-s/v w ^ = A* = 0.1-1 and e B = 
0.002-0.25 as typically inferred from broadband modeling of 
cosmological GRBs (Panaitescu & Kumar 2002; Chevalier, Li 
& Fransson 2003; Yost et al. 2003). Figure 3 shows that four 
of the SNe (19941, 1985F, 1997X, 1990B) have contours that 
independently rule out this entire parameter space while eleven 
place significant constraints. The late-time radio emission of 
SN 2001em is consistent with a typical GRB jet only if the mass 
loss rate is M- 5 /v w ^ w 1 (for e B « 0.2). Continued radio mon- 
itoring of SN 2001 em will reveal the origin of the emission, 
which may simply result from a late peaking radio SN within a 



dense circumstellar medium (Williams et al. 2002). 

5. DISCUSSION AND CONCLUSIONS 

Late-time observations of SN 1998bw (t ~5.6 yrs) have al- 
lowed us to test the hypothesis that GRB 980425 was a standard 
GRB viewed far away from the jet axis. Our measured upper 
limits at 1384 and 2368 MHz are consistent with the continued 
power-law decay of the SN emission. These limits imply an 
off-axis jet is only plausible if the normalized mass loss rate of 
the progenitor star is M-s/v w j < 0.04 (for e B > 0.1). This is 
~ 20-200 times smaller than the observed mass loss rates for 
local Wolf-Rayet stars (Cappa, Goss & van der Hucht 2003) 
and is below the range typically observed in GRBs. Larger 
mass loss rates are possible but only if the energy fraction in 
magnetic fields is low (i.e., e B <J 10~ 3 ). 

Even tighter constraints are derived for the off-axis jet model 
when we examine a larger sample of local Ibc SNe. The 
low luminosity limits derived for this sample require values of 
M-s/vw.3 « 0.01-0.1 or e B £ 10" 3 which are below values for 
typical GRBs. The absence of any late-time radio emission can 
therefore be used to put a limit on the fraction of core-collapse 
SNe that produce collimated, relativistic outflows. Our results 
imply that off-axis jets from nearby SNe are rare (^6%) with 
the possible exception that the radio emission from SN 2001 em 
is due to a GRB jet. 

This conclusion complements the findings of Berger et al. 
(2003a) who constrained the GRB/SN fraction through a radio 
survey of local Ibc SNe at early time. Berger et al. (2003a) used 
early, bright radio emission as a proxy for relativistic ejecta, as 
in the case for SN 1998bw. After studying 33 local SNe with 
detection limits 10 3 times fainter than SN 1998bw, Berger et 
al. (2003a) found no evidence for relativistic ejecta in any of 
the SNe observed, thereby constraining the GRB/SN fraction to 
£3%. 

Taken together, these results support a view that SN 1998bw 
was a rare and unusually energetic SN - distinct from lo- 
cal SNe and GRBs. In this scenario, the characteristics of 
SN1998bw/GRB 980425 are not dictated by the observer's 
viewing angle, but rather by the properties of its central en- 
gine. SN 1998bw was an engine-driven explosion (Li & Cheva- 
lier 1999), in which 99.5% of the kinetic energy (~ 10 50 erg) 
was coupled to mildly (T w 2) relativistic ejecta (Kulkarni et al. 
1998), while a mere 0.5% was detected in the ultra-relativistic 
(T « 100) flow. In contrast, GRBs couple most of their en- 
ergy to relativistic 7-rays. The observed diversity of cosmic ex- 
plosions (SNe, X-ray flashes, and GRBs) may therefore be ex- 
plained with a standard energy yield, but with a varying fraction 
of that energy given to relativistic ejecta (Berger et al. 2003b). 

We thank Edo Berger, Sarah Yost and Eli Waxman for help- 
ful discussions. AMS is supported by the NSFGRFR 
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Fig. 1. — The radio luminosity evolution of SN 1998bw is shown with late-time detection limits for 12 nearby Ibc supernovae 
(triangles) and three late-time detections for SN 19941, SN 2002ap and SN 2001 em. The four shaded regions represent an off-axis 
GRB jet for A* = M-s/v w ^ = (0.01,0.1, 1, 10) and assuming a range of values, e B = 0.002-0.25. We have assumed typical values 
of £51 = 1 and e e = 0.1. With the exception of SN2001em, all of these late-time observations are significantly fainter than the 
M_5 /v Wj 3 = 1 off-axis jet prediction and nine of the SNe constrain M_ 5 / v w ^ < 0. 1 . 
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Fig. 2. — Constraints on the off-axis GRB jet model for SN 1998bw derived from radio observations. Every flux measurement allows 
us to exclude the region of e B -M space rightward of the corresponding contour since this region produces a jet which is brighter 
than the observed radio emission at that epoch. The break in each contour marks the e B and M-s/v W 3 values that predict a peak jet 
luminosity equal to that observed. Above the break, a contour constrains a jet with rising flux, while below the break it constrains 
one with decaying emission. The sum of the excluded regions rules out the area rightward of the thick solid line. The region of 
reasonable parameter values for the putative jet is bracketed by M-s/v w ^ = 0.04-6 as derived from modeling of SN 1998bw radio 
light curves, and by e B = 0.002-0.25 as inferred for cosmological GRBs. Our observations of SN 1998bw allow us to rule out the 
majority of this region. 
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Fig. 3. — Constraints on the off-axis jet model for GRBs derived from late-time observations of local Ibc SNe. Our observations 
allow us to exclude the region of e^-M space rightward of each SN contour since this region produces a jet which is brighter than 
the observed limit. The region of reasonable parameter values for an off-axis jet is roughly bracketed by Ms / v w> 3 = A* = 0.1-1 and 
€b = 0.002-0.25 as typically inferred inferred from broadband modeling of cosmological GRBs. Four of the SNe contours rule out 
this entire shaded region while eleven place significant constraints. The recent detection of SN2001em is consistent with a typical 
GRB jet only if the mass loss rate is M_ 5 / v w> 3 « 1 (for e B w 0.2). The low end of the local distribution of observed Wolf-Rayet mass 
loss rates is also shown. 



